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Based on Childs finite length solution (ref. 1) for annular plain seals an ex- 
tension of the bulk flow theory is derived to calculate the rotordynamic coefficients 
and the leakage flow of seals with parallel grooves in the stator. Hirs turbulent 
lubricant equations are modified to account for the different friction factors in 
circumferential and axial direction. Furthermore an average groove depth is intro- 
duced to consider the additional circumferential flow in the grooves. Theoretical and 
experimental results are compared for the smooth constant clearance seal and the cor- 
responding seal with parallel grooves. Compared to the smooth seal the direct and 
cross-coupled stiffness coefficients as well as the direct damping coefficients are 
lower in the grooved seal configuration. Leakage is reduced by the grooving pattern. 


INTRODUCTION 


An important assumption for the reliability of high speed centrifugal pumps is 
a good rotordynamic behavior. Connected to this problem hydraulic forces acting on 
the rotor are of major importance. It is well known that neck or wear-ring seals as 
well as interstage seals (fig. 1) may have a large influrence on the bending vibra- 
tions of a pump rotor. Besides their designed function of reducing the leakage flow 
between the impeller outlet and inlet or two adjacent pump stages, respectively, the 
contactless seals have the potential to develop significant forces. This type of 
forces created by lateral rotor vibrations can be described by stiffness-, damping- 
and mass coefficients in a linearized model 
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For rotordynamic calculations of multistage pumps the machine designer needs to 
know this dynamic characteristics of the actual seal configuration. For smooth seals, 
where both stator and rotor elements have the same smooth surfaces, analytical and 
experimental investigations have been carried out (ref. 1 ,3,4,5) . The results con- 
firm the validity of equation (1) and dynamic coefficients can be predicted by the 
finite length solution derived in ref. 1 with sufficient accuracy. A finite length 
solution is also available for seals with different but directionally-homogeneous 
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surface roughness for the rotor and stator elements (ref. 2). 

The subject of this investigation is a seal type with parallel grooves in the 
stator element and a smooth surface in the rotor element. Concerning the leakage 
flow, this type of seal is more effective than the pure smooth seal, because of 
higher friction in the axial direction. However, if pump rotordynamics is important 
the dynamic seal coefficients have to be known. Up to now the stiffness and damping 
characteristics of grooved seals and their influence to the stability and unbalance 
response of pump rotors is not well investigated and there is a need for additional 
research in this area. 

At the present time, there are only a few techniques available in the technical 
literature. Black and Cochrane (ref. 6) have improved their earlier theory for smooth 
seals by introducing an equivalent length for the grooved section to reduce the cir- 
cumferential pressure gradients. Recently, Childs and Kim (ref. 7) have extended 
their analysis procedure (ref. 2) to predict rotordynamic coefficients of grooved 
turbulent annular seals. 

In the present paper an extension of the bulk flow theory (ref. 1) is given, 

to calculate rotordynamic coefficients and leakage for seals with parallel grooves 

in the stator. The theoretical results obtained by the developed procedure are cor- 

related to experimental results, measured at a seal test rig. Furthermore the grooved 
seal results are compared with corresponding data of the smooth seal configuration. 


BULK FLOW MODEL FOR SEALS WITH PARALLEL GROOVES 


Seal geometry 

Fig. 2 shows the type of seal, which is considered in our investigation. It 
consists of a smooth rotor and a circumferentially grooved stator. The seal has the 
radius R, the length L and may have different clearances at the entrance C and the 
exit C^ , respectively. The groove geometry is described by the groove deptR H R and 
the groove length L^ and land length L^. We assume, that the groove depth Hp has 
approximately the same order of magnitude as the seal clearances C , . 

In the following derivations the groove geometry is described simp?y by an average 
value Hp* for the groove depth (fig. 2). 


Bulk flow velocities 

In Childs finite length analysis (ref. 1) for plain seals a bulk flow model was 
used. Following this procedure we introduce the bulk flow velocities U„ in the axial 
Z-direction and Uq in the circumferential direction. The axial velocity U^ is con- 
sidered only in the region of the actual seal clearance H. Although there is a fluid 
circulation in the grooves: in the Z-direction, this part of the flow is neglected 
in our model (fig. 3). In circumferential direction the real velocity distribution 
(fig. 3) is replaced by a constant bulk flow Uq, which is assumed to act in the area 
of the average seal clearance H* = H + H *. At the rotor surface the fluid velocity 
is Ru) with the shaft angular velocity w . 

Fig. 4 points out the variables of the bulk flow model for a seal location with 
coordinates Z and 0: the two mentioned velocities U , Uq, the fluid pressure p and 

Zj 
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the local clearances H, H* = H + H R *. All quantities depend on the coordinates Z 
and 0 and the time t, as well. 


Wall shear stresses 


Concerning the wall shear stresses at the rotor and the stator we refer to Hirs 
formulation, expressing the shear stresses by an empirical function of the bulk flow 
velocity relative to the wall. First we apply this relation for the rotor assuming 
a smooth surface in the two directions. With the bulk flow velocity V R =[(U„-F.w ) 2 
+ u 2 ]l /2 relative to the rotor surface (fig. 5 ), we obtain R y 

Lx 


2H V r 


m 


= n T 


R P ? P P 

2 V R = C R 2 V R 


n R , m R are empirical turbulence coefficients, p is the fluid density and v the kine- 
matic viscosity of the fluid. 

With the relations of figure 5 we can determine the components of the rotor wall 
shear stress 


X R0 = T R (U 0 Rw)/V r 5 t rz = t r U z /V r 


(3) 


Contrary to the definition of t r „ in equation (3) we use a slightly different form 
and define x R with the seal clearance H* instead of H, which is used for t in the 
axial direction. After some steps we obtain from equations (2) and (3) RZ 
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with the axial Reynolds number R = 

a v 


Caused by the parallel grooves the stator has different friction characteristics 
in the two directions. In axial direction it behaves like a rough surface. For the 
circumferential direction we assume a smooth surface. Again we refer to Hirs formu- 
lation and express the wall shear stress t in dependence of the bulk flow velocity 

v s = (u e + u z )T 


- c It V 

' L s 2 V s 


(5) 


In the case of the grooved stator the triangle ratios of figure 5 are not quite 
correct. Nevertheless, we still use the relations and express the shear stress compo- 
nents approximately by 


S0 


= T S U 0 /V S 


( 6 ) 


T sz = T s V V s 
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The resultant friction factor Cg (eq. 5) can be introduced by the following 
superposition rule 


C„ = 


'SB 


22 221/2 
(Q^ cos cp + Cg Z sin cp) 
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= n 
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nig Z are empirical turbulence coefficients. 


Formula (7) describes the change of the resultant stator friction factor Cg in de- 
pendence of the local flow angle cp (figs. 5,6). If we consider one of the special 
cases, e.g. a pure flow in axial direction, we obtain 


^Z = ^S’ Uq = 0, cp = it/ 2, sincp = 1 , coscp = 0 
C S = C SZ i T S0 = °’ T SZ = C SZ T U Z ' 


From equations (5), (6), (7) we finally end up with the two shear force compo- 
nents 
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Cg', Cg are defined in the Appendix. (8a) and (8b) can be compared with prior re- 
sults from Childs and Kim (ref. 2). They differ only in the coefficients C ' , C ". 

O O 


BULK FLOW MOMENTUM AND CONTINUITY EQUATIONS 


Figures 7 and 8 show a differential element of the fluid having the dimensions 
Rd0, dZ, H(Z,0,t) or H*(Z,0,t) respectively. The upper and lower surfaces cor- 
respond to the rotor and stator seal elements and have the velocities Rm and zero • 
Figure 7 points out the bulk flow velocity components U y and U„ with their changes 
in axial and circumferential direction along the element. For the derivation of the 
momentum equations the wall shear stresses t , T and the pressure p at the different 
seal surfaces have to be taken in account (fig. §). Summing forces in the two di- 
rections for the free body diagram leads to the axial and circumferential momentum 
equations 9 and 10, respectively 
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In comparison to the derivations of reference 2 both equations have an additional 
term with the average groove depth H R *, expressing an added momentum change caused 
by the grooves. Furthermore the shear stresses are different as described in equa- 
tions (A) and (8). Note that H* is used in the circumferential momentum equation. 

The bulk flow continuity equation (11) also has 


3H 1 3 HU„ 3 HU 1 3U„ 

3t + R ( 30 1) + ( 3Z _) + R%* 30~ = 0 


an added term resulting from the flow difference in circumferential direction in 
the area of the grooves. If we substitute the shear stresses , t , t q, by 
the velocity-dependent formulas (A) and (8) we obtain the complete bulk flow equa- 
tions (see Appendix A), which can be used for further analysis. By introducing the 
following variables 
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the equations can be treated also in nondimensional form. 


PERTURBATION ANALYSIS 


In the further analysis we follow strictly ref. 2. The governing equation (9), 
(10), (11) or the corresponding equations in non-dimensional form define the bulk 
flow velocity components Uq, u and the pressure p as a function of the variables 
R0, z and the time t. The expansion of this equations in the perturbation variables 
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tion equations. 
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Zeroth - Order - Equations 


The eccentricity ratio £ = 0 describes the centered position of the rotor in 
the stator element. If we introduce u , u Q , h and p Q in the nondimensional bulk 
flow equations, we obtain the following zeroth order perturbation equations. 

Axial Momentum Equation: 


3p 1 

— — = - — — |o„ a „ + a a De + Aq} 
~ 0 „3 R oR sz osz 

3z 2 f 


Circumferential Momentum Equation: 


(13a) 


3 Up. 1 h * m h * m 

3^= " 2f ^ ( TT _) a R a oR (u 0o _1) + ( "h } " °s0 a os0 Da u 0o^ 


(13b) 


Contuniity Equation: 



The parameters of these equations are defined in Appendix B. h = f is the dimen- 
sionless clearance function for the centered position and q = ?C ) / (C +C* ) is a 
measure of the degree of taper in a seal. For the constant clearance seal, treated 
later in this paper, it follows q - 0 and h Q = 1 . The quantities a R , o^ z and °S0 
are defined by 


°R = ‘C 1 X R ; °sz ' ( C 5 X sz ; °s0 = ( C } X s0 
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The solution of the zeroth order equations define the steady state leakage and 
the development of the circumferential velocity u (z) due to wall shear. In general 
the coupled and nonlinear equations have to be solved iteratively to determine the 
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leakage flow corresponding to a specified pressure drop. A leakage coefficient C 
can be introduced for the leakage/pressure drop relationship 


A P„ 



2 


In formula (16a) the pressure drop at the entrance 


(16a) 


(,6b) 

(1+q) 2 

is included. For the special case of a constant clearance seal without fluid rota- 
tion we obtain the simple relation 
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3 z 


( °R 


+ o 


) 

sz 



(17) 


which can be used to determine the empirical coefficients n^, m„, n , m 

R’ R’ sz’ sz 


First order equations 

The first order equations describe the pressure and flow conditions due to a 
small seal motion about the centered position. Their derivation is relatively ex- 
tensive, therefore we mention only some important steps in the solution procedure, 
following again references 1,2. To find results for the first order quantities u . , 
U Q1 , P-) > the time and 0- dependency is eliminated by the assumption of a harmonic 2 
pressure and velocity distribution in circumferential direction and by introducing 
a. circular harmonic seal motion with the relative radius r = R /C and the frequen- 
cy fi. In this way the first order equations are reduced t8 a system of three cou- 
pled, complex ordinary differential equations for the complex unknowns u , Uq and 
p.| , which now depend only on the axial coordinate Z. 
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The 3x3-matrix A has the following elements 


(18) 
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with r = ft-o)U Qo (z), T = L/V 
and the right hand sight consists of 
§1 = 2q/f 3 + jrT/f 

S 2 = - f * A l0 

g 3 = -(A 1z + 2q/f 4 + jrT/f 2 ) 

The parameters A 10 , A^g , A^g , A 1z> A 2z , A 3z are expressed in Appendix B. 

With the boundary conditions of ref. 1 the equations (18) can be solved and 
yield the solution for the velocity and pressure field of the form 
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Dynamic Coefficients 


From the pressure field solution of (19) the reaction force components acting 
on the rotor due to the circular shaft motion can be determined by integration of 
the pressure along the seal and in circumferential direction. As pointed out in ref. 
1 the nondimensional form of eq. (1) 


J 

nRAp 



( 20 ) 


can be used in the definition for the radial and circumferential components of the 
reaction force 
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Finally the dynamic seal coefficients K, k, C, c, M, m can be found by a least 
square curve-fit-procedure applied to the right hand side of equation (21). 


Applicability of the derived equations 

The presented equations can be applied for smooth as well as grooved seals. 
Furthermore it is possible to investigate constant clearance and tapered seals. All 
equations are expressed in a form corresponding to prior derivation from Childs 
(ref. 1,2). Therefore it is easy to reduce the general form into simpler expressions 
and to compare them with prior equations. In the special case of a smooth constant 
clearance seal it follows H„ = H* =0;H=H*;n o =n =n ; m„ = m = m : 

C = C" = 1, q =0 and h = R 1 . B B 32 » R 32 33 

All further parameters, especially those defined in the Appendix are changing for 
this special case and the resultant equations coincide with the equations in ref. 1 . 


THEORETICAL AND EXPERIMENTAL RESULTS 
FOR SMOOTH AND GROOVED SEALS 


The objective of the theoretical and experimental investigations is to check 
the usefulness of the developed model for grooved seals. For this task predicted 
and measured dynamic coefficients for grooved seals are compared. It is of further 
interest to point out differences of leakage and rotordynamic coefficients for smooth 
and grooved seals. 


Geometry of the Test-Seals 

The two seals which have been investigated are shown in fig. 9. Seal 0 is a 
smooth constant clearance seal without any grooves, seal (2) has eight parallel 
grooves in the stator. The rotor elements are considered to be smooth. Both seals 
have the following data: radius R = 23,5 mm, length L = 23,5 mm, constant clearance 
C = 0,2 mm. For the grooved seal the groove depth is H = 0,5 mm with L r = 0,7 mm 
and L = 1 ,5 mm. M ^ 


Test Rig for Seal Investigations 

With the test rig shown in fig. 10 the leakage flow as well as the dynamic seal 
coefficients can be determined by measurements. The cross section shows two seal 
inserts integrated symmetrically in a very rigid housing. A stiff shaft, driven by 
an ac-motor, rotates inside the housing and acts as the second part of the seal. 

The housing is flexibly supported by eight beamlike springs and therefore a pure 
lateral jnotion relative to the shaft is specified. In the operating condition water 
with 30 C is entering the housing in the center, flows througp the two test seals and 
is exiting the housing at both ends. 

To characterize the fluid state several pressure and temperature pickups are 
distributed at the test apparatus. The fluid velocity, determined from the mass flow 
rate, was measured in the supply line. The housing can be excited by test forces, 
which are measured by a force transducer. To detect the resultant motion of the 
housing relative to the shaft eddy current-pickups were used. 
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For the steady-state leakage measurement the housing is fixed in a centered 
position and the mass flow rate and the pressure drop are taken by measurements. 

With this test data and additional informations concerning the fluid the leakage flow 
and the empirical turbulence coefficients can be found. 

The experimental determination of the seal dynamic coefficients works with a 
parameter identification procedure. In the measurement step test forces are applied 
to the housing in radial direction and with the measured relative motions between 
the two seal surfaces mobility frequency response functions can be calculated by 
Fast Fourier Transformation. Corresponding frequency response curves of a seal test 
rig-model are fitted to the measured functions by variation of the seal dynamic coef- 
ficients. The identification procedure is described in more detail in ref. 4,5. 


Leakage Performance 


To compare the leakage flow of different seal stators the leakage coefficient 
C r is defined as follows 

Lj 


Q 


O / 2A r, 

C. 2ttR 2 / — 2-° 


( 22 ) 


Q is the volumetric steady state flow rate measured at the test rig and Ap the cor- 
responding pressure drop. C is a nondimensional relative measure of the leakage 
expected trough seals with the same radius. Fig. 11 illustrates measured values C R 
for the two seals without grooves and with grooves. The grooved seal (2) has a 
leakage coefficient, which is about 15 % lower than that of seal (]) without grooves. 
Both leakage coefficients increase slightly with the pressure drop Ap. 


Empirical Turbulence Coefficients 

In the described analysis the friction factors X R > X > X ~ were characterized 
in terms of the empirical turbulence coefficients (eq. 15f? These empirical coef- 
ficients have to be determined from static test data before a theoretical prediction 
for the seal dynamic coefficients is possible. Leakage rates and pressure gradients 
are measured for this task. The steady state axial pressure gradient is described 
by eq. (17) 


3P o P -2 

9 z~ ( °R + °sz ) 2 V 


(17) 


With the measured pressure gradient, the velocity V and the density p the combined 
0 -values in paran theses can be calculated. We start with the smooth seal configu- 
ration and suppose that o„ can be applied for both the smooth rotor and the smooth 
stator. From eq. (17) we obtain 2*a R and the friction factor X R , respectively. The 
second test is carried out with the smooth rotor and the grooved stator. From the 
measured quantities for this case first of all only the combined friction factor 
(d R + o ) is known. With the value a R from the first test a ^ and the corresponding 
X^ z can be calculated. Fig. 12 shows he two friction facto? z X R and X independence 
the axial Reynolds number R . The grooved seal (g) has approximately tne double 
friction factor compared to a ?he smooth stator of seal © • 


138 


From the \ versus R data the empirical turbulence coefficients n , m R , n sz , 
m are calculated with a- least square procedure, based on equations (T5). 

S7 The following values were obtained (fig. 14): 


n R = 0,062 ; m R = - 0,22 

n = 0,058 ; m = - 0,13 

sz sz 


(23) 


Concerning the values n , m , we assume, that the smooth surface constants are 
relevant. ^ ^ 


Dynamic Coefficients, Influence of the Groove depth 

With the presented analysis we can now calculate the dynamic seal coefficients. 
The numerical procedure is applied to determine especially the stiffness and damping 
coefficients K, k, C, c for seal (2) in dependence of the rotational speed and the 
average groove depth. The axial Reynolds number R is kept constant in this in- 
vestigation. Besides the seal geometry (fig. 9) further input data are as follows: 

Fluid data for water with 30°C 

Axial average velocity V =16,46 m/sec 

Entry swirl U = 0,2 Rco 

Inlet pressure loss ^ = 0,5 

Describing the friction behavior, the empirical coefficients from (23) are used 
(fig. 14). 

To point out the influence of the groove depth, different values of H R * are 
assumed (figs. 13,14). Taking the average depth from a geometrical approximation 
(sum of the upper areas equal the sum of the lower areas) we obtain H * = 0,1 mm. 

The selected values for the calculations are H R * = 0; 0,1; 0,2 mm. Note that the 
empirical constants are held constant in this investigation. 

Fig. 13 illustrates that the direct stiffness and damping as well as the cross 
coupled stiffness are reduced by H *, there is a weak influence to the cross coupled 
damping. An increase of H R * from 0 "to 0,2 mm reduces K about 50 % and C about 37 %, 
respectively. 

In Fig. 14 the direct stiffness K and the direct damping C are compared for the 
two seal types (]) and (g) . All presented values are related to the coefficients K 
and C of seal ( 7 ) without grooves. The results correspond to R = 8419 and a rota- 
tional speed 4000 rpm. All other data are the same as in the example before. Starting 
with the reference seal (T) with empirical constants for smooth surfaces the stiff- 
ness ratio as well as the damping ratio are equal to 1 by definition. Seal (2) with 
grooves has other constants (fig. 14). If we consider only this influence (change in 
friction) and keep the average groove depth constant (H R * = 0), direct stiffness as 
well as damping increase to 1,12 and 1,19 respectively. Taking now the empirical 
values of seal (g) constant, an increase of H R * reduces both stiffness and damping, 
as was shown already in fig. 13. If we compare with corresponding measurement re- 
sults, we recognize that the calculation with H R * =0,1 mm yields good results for 
stiffness but 15 % too high values for the damping. 
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Dynamic Coefficients, Comparison of Test Results 
and Theoretical Results 


As described before the dynamic seal coefficients can be measured at the seal 
test rig. For one working condition with constant rotational speed, axial velocity 
and constant fluid temperature four frequency response functions are measured by 
exciting the dynamic seal test rig and picking up the response in the two directions. 
A computer takes over the measured data and calculates the dynamic seal coefficients 
by means of a least square identification procedure (references 4,5). Several meas- 
urements were carried out for different rotational speeds but constant axial Reynolds 
number R and fluid temperature. Fig. 15 illustrates for the two seals with and 
without frooves the identified stiffness and damping coefficients versus the rota- 
tional speed. The test results show, that neither the direct coefficients nor the 
cross coupled coefficients are equal in amount, as expected from theory (eq. 1). 

Each of the two coefficients, which should be equal in magnitude, are shown in the 
diagram as found out by the identification process. The scatter of the measurement 
points is very different. The added mass terms are much higher than predicted by 
theory. They are not presented in the diagram. 


The fitting curves show the tendency as expected from theory, a slightly para- 
bolic decrease of K with the rotational speed, constant direct damping and a linear 
increase with the rotational speed for the cross coupled terms k and c. 


Seal (7) without grooves has higher direct stiffness and damping terms compared 
to seal (D with grooves. The cross coupled stiffness of seal (T) is also greater 
than of seal (g) . There is no clear difference in the cross coupled damping values. 
Basicly, the different measurement results for the two seal types show the expected 
influence of the grooves (see also fig. 13), reducing especially K, C and k. 


In fig. 16 the identified stiffness and damping coefficients of seal (2) are 
refered to the corresponding predictions from the grooved seal model. There are two 
parameters in the model, which can be changed slightly, to obtain a better correla- 
tion between measured and theoretical results: the average groove depth Hp* and fluid 
entry swirl at the seal entrance. It was found, that values for the entry swirl be- 
tween U„ =0,1 Rw for low rotational speeds and Up = 0,3 Rw for higher rotational 
speeds and a average groove depth of H R * = 0,1 mm were best suited to fit the theo- 
retical values to the test results. The predicted direct stiffnesses are slightly 
lower, the direct damping values about 20 % higher than the measured quantities. 


The influence of the axial Reynolds number for a constant rotational speed of 
4500 rpm is pointed out in fig. 17. Measured as well as predicted parameters show 
the similar trend for the different seal coefficients. K is increasing in a para- 
bolic curve, C depends linear on the Reynolds number and the cross coupled terms k 
and c behave indifferent. 


Finally in fig. 18 the seal model predictions for the two seal types are com- 
pared, in accordance to the measurement results in fig. 15- The results are as fol- 
lows, an expected decrease of both stiffness coefficients from seal (7) to seal (2) 
but a weak difference of the damping parameters for the two seals. 
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CONCLUDING REMARKS 


A theoretical model, based on Childs finite length solution, is presented to 
determine rotordynamic coefficients and the leakage flow of seals with parallel 
grooves in the stator. Calculated and measured stiffness and damping values for the 
investigated seal with eight grooves show, that the developed model is useful for 
the prediction of this seal type. 

Concerning the pump efficiency seals with grooves have the advantage of a lower 
leakage flow compared to the smooth seal. But the direct and cross coupled stiffness 
coefficients as well as the direct damping are reduced in the grooved seal configu- 
ration. This has to be considered when rotordynamic problems are important. 
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APPENDIX A: COMPLETE BULK FLOW EQUATIONS 
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APPENDIX B: PERTURBATION COEFFICIENTS 
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The parameters D g , D a , D ei ,D ai are very extensive and therefore not presented. 
If desired, they can be obtained by the authors. 
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higure 1. - Seal types in turbopumps 
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Figure 2. - Seal with parallel grooves in the stator 
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Figure 3. - Bulk flow velocities 
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AXIAL DIRECTION 


Figure 4. - Variables of the bulk flow model 




Figure 5. - Wall shear stresses 
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Figure 6. - Stator friction factor 
in dependence of the flow angle 
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Figure 7. - Bulk flow velocities at differential element 
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Figure 8. - Control volume for a seal element with 
pressure and shear stresses 
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Figure 9. - Investigated seals with and without grooves 
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Figure 11. - Leakage coefficients for seal 1 
and seal 2 
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Figure 12. - Friction factors for the stators of seal 1 and 
seal 2 
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Figure 13. - Calculated stiffness and damping coefficients of a grooved seal 
in dependence of the rotational speed and the average groove 
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Figure 14. - Comparison of calculated and measured direct stiffness and 
damping coefficients for seals 1 and 2 



Figure 15. - Measured stiffness and damping coefficients in dependence of 
rotational speed for seals 1 and 2 
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Figure 16. - Measured and predicted stiffness and damping coefficients in 
dependence of rotational speed for the grooved seal 2 



Figure 17. - Measured and predicted stiffness and damping coefficients in 
dependence of the axial Reynolds number for the grooved seal 2 
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Figure 18. - Predicted stiffness and damping coefficients in dependence 
of rotational speed for seals 1 and 2 



